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VIII. — Researches in Physical Astronomy. By J. W. Lubbock^^ Esq. V. P. 

and Treas. R.S. 



Read February 9, 1832. 

AN general^ two methods present themselves of solving any mechanical pro- 
blem : the one furnished by the variation of parameters or constants, which 
complete the integral obtained by the first approximation; the other furnished 
by the integration of the differential equations by means of indeterminate 
coefficients, or some equivalent method. Each of these methods may be 
applied to the theory of the perturbations of the heavenly bodies ; and they lead 
to expressions which are, of course, substantially identical, but which do not 
appear in the same shape except after certain transformations. 

My object in the following pages is to effect these transformations, by which 

their identity is established, making use of the developments of R and r -j^ 

given in the Philosophical Transactions for 1831, p. 295. The identity of the 
results obtained by either method serves to confirm the exactness of those 
expressions. 

Integrating the equation 



d^r^ _±_ ^ [^ 



+ 2/dE + .(^)=0 



2dt^ r a 

omitting the terms which are independent of the quantities &, and which 
result from the part of R which is equal to ^ ^^3 -^, and the factor ^. 



a 






— _ I Z>, . -f- —-i: > cos i (w e -— ni) 

^ IfJ- Wl } ^^''' C ^"^ ^ "^ ''' + w f- tor) 
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n^ /i^A _ (1 4-21) gg. 



i (« ■— n^ + 2 W^ I (« — «^) 



a 3 ^»* 



-y- 



3 {i(w — w^) + w}2y2 J n a. 

■ i n .. — Tj . r? .. .j.i iii i i .w ii rw.rH iW wr.i M WWM H" ii ■ " ■ i J — . . . i m -n- , ■_ . . . .■ i .i . .i . i- — i.,iii «„ i „ .,.. », , „ .. „ ..i.- „ . , ,^_ . , , „ j.„ , .. .^ ^. w ^i — > i.ii r ii »r .iii.inii« i | Pi In i i Mi. ai ii « MW ffi l.i ■ r . f . mm S^ «aiw - i .. w i » — — »ii.- —.... — fj . t 

2(1 (n - «^) + 2 ?a) i (w - n,)(i (n - ?2^) - 91) ^ ^^ ~ ^') ^^ "' 

2a^daJ \^ '^ / 

^i (^i — n^) 4- f2^^ M (w — «^) + ?i 4- w^^ f i (w — w^) — ?a + w^ ) ^' ' ' 

, ^^ /3 (1 + i) «% _ i« . 

i being any whole number positive or negative^ but excepting the arguments 

n^ ^ n J n n 1 

*~»««««™--»««*w--«**»««*»>*»*»-^^ BMMM 'i w ii- ww i*... — i^. i I — ^ ■ L I i n ii . i n iwi i iii r n ■j iiii w iiiiiii ii i HHUmn — mm ii w ■w u i ■i i iU MJ I H ii m i '. iilii Hi ii M* wl innKi - '■■•a* " " —" ■ '"■ " i "" "" m i-iiiiWW » i iyi i ri i i ii m l mi t i i mi. ii jbl i , n i >■ w . wmnHfusu i , i»- ^ 

(« - ?i^) (i {n - w^) + w) (i (w - n^) - n) ^ - ^i L2 (i (?J - n^) - w) 2 (i (n -r- ?2^) 4- ^) -^ 



4 



2 (i (w - n,) - n) ^ ~" ^^ 2 (i (w - w^) 4 n 

ResolvinE' the other fractions in the same manner. 



r 



i^^- — b, . cos lint — nd) 



— — — — — — — -:-" s -^ ^1 ,• "- — . bo s 4- TT-—: ^ai 1 4- 7; — ha -11 ? COS % (n t — • n,t) 

+ ,^ ^__J — I i — ^_^ 5 — ^_^ J ^ .^ — . 5 V ecos ( z (?i i — Tzi) 4. n ^ ~~ tJT 

{i (;i - n.) 4 w} I 4 a^^ s»*-i 2 a/ ^>* 4 a^^ ^*'+^ J V 






.„^ J f 2 1 4 2) / ^^ jy 4, ^^ ^ _ 3a^ , 1 _ jg^ , 

2{i{n-^n^) + 2n}\^ l4a^2 3,i-l + 2a/ ^'^ 40^^ 3,i+i / 4^i%i-i 

, (I 4-2i) a3 Sia% 1 /. , , a , . \ 
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■ a»iyw> 



n e 



2i{n -- n^) 



^^ ' + ^^ {4^^ ^3,i-I +2^3 ^3,^ - IJ *3,i+l} + ifl \i 



(1 +2i) o^ , 3ia^ , 



+ 



i 






a ^Ki 



3n 



4i (n — n^) 



a. a^ da 






n 



iah,^ a ^h,. 



A(i{n-'n^-\-2nj^ ^' 



a^ da 



>e co%li{nt— n^t) + nt-^'m 



2in 



Q {n — n^) + n^J 



J3a^ ^ 

14 a;^ 3,8-1 



a 



5... 



2 a ^>^ 4 



^2 ^3,^+l } ^i COS (i (^ ^ - w^O + n,t - ^, 






n 



2 n (?2 — w^) — ?e + n^j 



7)H 



a 






i.i + 



2a^2 3,8+1 / 



3( 



n 






.e^cosn(w^ — w^^) + n^t—'ZAf^ 



2(i{7i--n^ "-n + n^j 



2i j 



3a 



a 



4^.^3,^-1 + 2^^3,8+4^^ 






^ 4 a^^^S,i-l ^^ %• —4— «;.%e+lj ' \ ^ 

serving that 

^'* 2a^ I ^'*~^ 3,84-1 J 



w^^) + w^^ — "zsr^ 



cfd5, .. 



a a^ L a^ ^»^ 2 



_ 1 , 

3,8-1 2 ^ 



,8+1 j 



the preceding expression may be put in the form 



a 
r 



r — r h, . COS i (w ^ — n.t) 



n 



2{i (n — n^) + n} 



1 2 a^2 ^3,8-1 ^^3 ^3,8 ^ 2 a^2 3,8+1 J '^ ^ ^ 



+ 



/ \\ 4a2 3,8-1 2a 3 2,8^4^2 3,8+1/ \^ ^ i/ 
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ne 



2{i(n—n^) + 2fi^ 
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f (i+2) a^ , (3 + 4i) a'^ r ___ (9 i + 6) a'^^ , 



cos 



(i(nt-- n^t) + w 1^ — -nr I 



+ ^, ^^ I (3i + 2) -^ b^. , + -^ 5., . 



(Si +6) a^ 



57 3,i+l/ 



cos 



n (w ^ — w^^) 4- w ^ — -nr I 






n 



2i(«-w^)^a, ^'^ (i(n^n,)-n') 



\257^ 3,i- 1 ^3 ^3,f 25;§ %i+ 



4i {n —n 



{2a^ y aP , a^ , \ 



iim| — 



n 



4(i(n-w^) +2n)^^' ^'' 



°^^ ^^a,- — ^ ^q,-_Li r Ucosf i(wJ — wi) 4- n^ — '^^■ 



ne, r(3 4-9i)«^ A 

(\ I 4 /? 2 3,2—1 



2ia 



^3,e —4 a/ 3»^*+l 



cos 



f i (w i — • n^^) + Uft '^'mA 



ne. 



2 (^H^ — ^/) — w + w^^ 



{ 



(3i-3) 



a/ 3,i-l 4 a,^ 3,2+1/ 



cos 



H (w ^ — n^t) 



+ w^^ — tzr^ 



) 



and by further reductions 



a 

r 



nab, . 

1)2 






a" 7 
" 3 S,J 



+ 



lr^^3'^+i}} 



cosi (wi^ — n^t) 



+J 



(i+l)« 



A (n — Til) + nj 



I 4 a/ 3>*- 






1 2 a/ 3»* • 4 



+ 



n 



(i(w — w^) + 2wy 



L 8 a/ 3''-l ^ '^ a,3 3,1 4 a^s 3,,+!/ 
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Sw^a ^i,f 



2i(?i — w^)^ a 



S-ecoB 



li{nt-^n^t) -{- nt — 'm\ 



+ ^ 



n 



(i(„_.,)+„,)l4«.^ 3,i-l 2 



f3a2 , a , «2 , 1 



w / (3 4- 9 i) a^ ^ 

(i (n ^ n^) + n + n^j ^ 



^a-S,?--! J-^Sd 



ta 



(I -^ i) a^ J 1 

"~8~~ ^7^ ^^^"+1 J 



n 



(j (n - K^) 



/ (3i-3) ga _(3t-l)a^, "1 



e^cos 



t being, as before explained^ any whole number positive or negative^ excluding 
only certain arguments, 0^ nt + B — isT^ and nt + b--^^. 
Considering the terms which have hitherto been neglected^ if we suppose 

-^ = 1 + Vq + ecos (n (\ + k) t + E —w\ + ej^ cos ln{l ■{- h^) t + s --^ 'm\ 



we have 



a 



""^ - 2^ ^^'^ 



02 



2a^2 3,1, 









See Phil. Trans. 1831^ p. 53. 
If ?z (1 + 2 r^) = n and n^ = -^ if e is the coefficient of sin (/i # + g — m) in 
the expression for the longitude, and/] is determined so that the coefficient of 
sin {nt '\' B — '^) in that expression equals zero. 



T='-4'"+Tfi7.*- + 4'+4--'"«-5l^''")'»'(»('-^A.)'+ 



S -— -37 



+ e^ 



{1$ ^^'^ - T^, ^'•' + 8^^^'4 '' '"' (» (1 + *-) * + ^ - --) 



m. 



In the theory of the moon replacing ~~^, 



a 



r \2i^a3 



' + e 






3 m, g3 

4 /uc- a^3 



J i^ + £ — - 'ZtT I 
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The preceding results are obtained by the direct integration of the differeii« 
tial equations : I shall now show that they coincide with the results obtained 
by the variation of the elliptic constants. 

The equations for determining the variations of the elliptic constants are^ 

d jR 

d a = — 2 a^ w -r~- d ^ 

as 

11 e— \ r ^ ae da 



an i/\ — e^ ^^ .-----— .dE,, , an ^/l — e'^ dR . . 

Z. (I — Vi _ e2) -~j™ d ^ H — T— d t 



dg e d t3* 



, aw//l — e^d JJ , 

dtrss: «- T~ a t 

e de 

sini \/l — e^ d* 

an d jR 

■""sini ^/l — e^ dv 

See the Th^or. Anal. vol. 1. p. 3?0^ or The Mechanism of the Heavens., p. 23 L 
In these works R is used with a contrary sign to its acceptation in the 

Mecanique Celeste, which I have followed. 
When the square of the eccentricity is neglected in the value of the radius 

vector, the equations may be employed in the following shape : 

da:=^2a^^.4^d^ da~.d^=^^d^+2a^^^df d e = ^L!if 4^ d^ + - i? cU 

dfi ede da 2dg ed'zsr 



If ?==y\d^ d^ — ^jyandRdt 



a ^ a 

r a 



< 1 + e cos (w i + £ — -zo") > + cos (nt + S'-^m) Se — enin (nt-^ s — m) {$ s — $'m) 



4- 2 e cos (2 n ^ + 2 s ~ 2 ^) ^ e -- 2 e^ sin (2 7^ ^ 4- 2 s •— 2 ts*) (^ g — ^ ts-) — e sin (» e 4- s --. w) ^ 



a na 

r 



-------- b-, . cos lint --- n.t) 



< cos (w ^ — -sr) cos l^w^—- w^^) -j-w^— ts-l 
+ sin (nit — -sr) sin ({{nt-^-n^t) 4 w ^) — ikt | > 



(i-{- l)fi 
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/_^^^. , 4- _^5,.-~|-^&<,. , 1 \ cos (iint-^-n J) +nt'^'m) 



2ne / _ (2 + ») g^ 7. _ iL±Jl SLb +(?L±lIl^b \ 

I COS (n t — 'm) C08 li (n t — n^ t) '^2nt-^ 2m\ 
+ sm(nt -^ 'm)sm H (« ^ — w^ ^ +• 2 w ^ — 2 ra- j > 






we f 3a^ , _ "^ 5 _ J?L 5 l 

{i(„_«^)_«} l4a,« 3.^-1 2V ^'' 4a,« 3.'+M 

< COS (2nt-'2w) cos n(»^— w^^) — w^-fisrj 
— sin {2 n t •— 2 'm) sin ( i {n t -^ Uf t) — n t + m\ 

nae ^1,£ ^ ^ mna o^, 



4 (?z — ?ey) a^ 



y \ w /7 p 1 i / \ 

cos I i (» ^ — w. + w ^ -- ^ I — ^ — - — '- cos I i (« ^ — n, e) 4- ?i ^ «- -nr I 



a^ ?J g ^ - '^ "^ ^ vi^ n o ^1 »• 



(?? — ?2,) a 



_^ii!. cmli(nt'--nj) + nt^w) -— ^ . , .- -~^ cos { i(n t -- n^t) -j- n t -^ t^ ] 

da \ ^ ^ ^ / 2z(n— w^)^ «i \ / 



{ i (^ _ n,) + ?? + 71 J. L "~8 a;^ S,e - 1 ^^ 3,e 3 ^^2 S,e+ 1 / 

< COS (w i5 •— -^o*) COS n (w ^ -- ?2^ i5) + w i --- -sr + f?^ i — -sr A 

— sin (» i5 — -kt) sin I i (w ^ — ?z^ 1^) + w ^ — to* + ?i^ ^ — -sr^ | > 

{i(^_^^) _^4.^J 1 8 a/ ^'^-^ 8 a/^3,i4-i/ 

< cos (nt — 'm) cos H (« i5 — «^ 1^) — « 1^ — -zr + ;i^ i — • -sr^ | 

— sin (w 1^ — -sr) sin H (m ^ — w^ i) •— w ^ —• -sr + 71^ 1^ — sr^ | > 

It is easily seen that this expression is identical with that of p. 2325 obtained 
by the direct integration of the differential equation of the second order. 

Considering the arguments 0, ?2 ^ + ^ — '^5 w i? + 2 -- ^p stilly however^ neg- 
lecting for an instant the term ~i b^^^ e e^ cos {m — ^^) which requires parti« 

cular attention^ 

2 H 2 



a 
r 
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_ = ~ — - 63 A — -Ti — ;; Oo 1 + < • — O3 — -i — ^ 0? 1 > e COS (?^ ^ + g — tS*) 

4 a;^ ^'V / 

The terrti ^^ &g .^ ^ ^i cos (ts* — sr^) in the development of R gives 

d?<T= — -J — 63 06, cos (w — -(xr,) de = — - — - h. oe.sin te — w,) 

If /irsesinzsr Zrsecostar /i^= e^ sinter^ l^=: e^co^ 'ui^ 

d/i = ecoswdt3*+ smwde = — -— - 63 oe,cos w,= — ^i-^ 5. ol, 



d^ss — esinizrdOT + cosOTde = — ^^ 63^^ e^ sin ts*^ = — ^~^^ ^3^^ ^^ 



4a^2 '*^ ' ' 4a/^ 

The integrals of which equations are 

h = iVsin (gt + C) 1= Ncos (gt + C) 

h^=zN^sin(gt + C) l = N^coB(gt + C) 

Ng^^^J,,,N, isr,g=— J5i53^^IV 

e cos (w < + s — -sr) = iVcos (w — g) ^ + e — C) 

which will agree with the previous solution^ p. 233, if 

This theory of the sec alar inequalities appears to require to be extended to 
the terms depending on higher powers of the eccentricities ; but I may remark 
that the coefficient of the term e^ ef cos (2 tsr —- 2 w^ in the development of R 
vanishes in the theory of the moon^ or at least such part of it as is multiplied 



